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We present measurements of partial branching fractions of B — > K^ X, B — > K^X, and B — !> tt^X, 
where X denotes any accessible final state above the endpoint for B decays to charmed mesons, 
specifically for momenta of the candidate hadron greater than 2.34 (2.36) GeV for kaons (pions) in 
the B rest frame. These measurements are sensitive to potential new-physics particles which could 
enter the b — >■ s{d) loop transitions. The analysis is performed on a data sample consisting of 383 
xlO^ BB pairs collected with the BABAR detector at the PEP-II e^e~ asymmetric energy collider. 
Our results are in agreement with standard model predictions and exclude large enhancements of 
the inclusive branching fraction due to sources of new physics. 

PACS numbers: 13.25.Hw, 12.15.Ji, ll.30.Er 



B mesons decay predominantly to charmed mesons 
through the tree level process b c, while the tree am- 
plitude b ^ u and the one-loop processes b ^ s and 
b ^ d are strongly suppressed. In the standard model 
(SM), the inclusive branching fraction of B mesons to 
charmless final states is of the order of 2% [l[ . Particles 
associated with physics beyond the SM, such as super- 
symmetric partners of SM particles, could enter the loop 
amplitudes while leaving the tree-level processes nearly 



unaffected, making a sizable enhancement of the inclu- 
sive 6 — >■ s {d) g (where g denotes a gluon) branching 
fraction possible 0, S] • Additionally, since semi-inclusive 
processes are usually affected by smaller hadronic uncer- 
tainties than those that arise in calculations for exclusive 
final states, these decays can be sensitive to nonpertur- 
bative amplitudes, such as charming penguins 

An interesting theoretical mechanism that can mod- 
ify the SM prediction is provided by the Randall- 
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Sundrum framework, in particular from the Warped Top- 
Condensation Model where a radion field cj) is postulated. 
In the case where 1 < m{(j)) < 3.7 GeV, the radion would 
decay dominantly to gluons, thus enhancing the rate of 
the charmless B decays through the process b — >■ scf). In 
such a model the 6 — > s inclusive decay rate could be 
enhanced by an order of magnitude with respect to the 
SM predictions ^]. 

Historically, an enhancement of charmless B decays 
had been postulated \d\ to explain the deficit of 6 — >■ c pro- 
cesses observed by the ARGUS and CLEO experiments 
0. Later measurements and refined theoretical calcu- 
lations established that no significant discrepancy was 
present Q. Inclusive b ^ sg decays have been searched 
for by the ARGUS, CLEO, and DELPHI collaborations 
0. None of these experiments has found a statistically 
significant signal and only upper limits in agreement with 
theoretical expectations were set. 

In this paper we present measurements of partial 
branching fractions of inclusive charmless i?-meson de- 
cays. The signature of these decays is the presence of a 
light meson (if"*", if^, or 7r+ [13]) with momentum be- 
yond the kinematic endpoint for B decays to charmed 
mesons, measured recoiling against a fully reconstructed 
B meson. It is possible to compare our results with the 
inclusive branching fraction of 6 — )■ S7 in the same kine- 
matical region and with some recent theoretical predic- 
tions Q based on Soft Collinear Effective Theory. 

The measurement is performed on a data sample col- 
lected by the BaBAR detector [ll| , operated at the asym- 



tral pions, and at most two K'^ tt+tt . 
struct D*- Wit-; D*" D"n°; W 
K+n^Tr°, A'+7r^7r"7r+. K2n+n^\ and - 



We recon- 

^ K+TT-, 



metric energy e+e" PEP-II collider at the SLAG Na- 
tional Accelerator Laboratory. We use 347 fb~^ (equiva- 
lent to 383 X 10^ BB pairs) collected at a center-of-mass 
energy corresponding to the mass of the T(4S') res- 
onance, which predominantly decays to charged or neu- 
tral BB pairs; a smaller sample (37 fb^^) of data col- 
lected at an energy of 40 MeV below the T{AS) peak is 
used to study the background originating from contin- 
uum e+e^ ^ <N (q — u, d, s, c) processes. 

In order to achieve the highest possible control over the 
backgrounds (mostly arising from continuum), we fully 
reconstruct one of the two B mesons (denoted by i?rcco) 
and search for a high momentum light hadron [K^ , K^, 
or 7r+) among the decay products of the other B (Bsig)- 
The full reconstruction of the Breco candidate allows us to 
determine the four-momentum of i?sig precisely. In order 
to suppress backgrounds arising from the dominant B 
decays to charmed mesons, we require the light meson's 
momentum p* in the Bgig rest frame to be greater than 
2.34 (2.36) GeV in the kaon (pion) case. The separation 
of from 7r+ candidates is based on the Cherenkov 
angle measured in the Detector of Internally Reflected 
Cherenkov light. 

The Brcco is reconstructed in the decays B — > D'^*''Y^, 
where is a combination of hadrons containing one, 
three, or five charged kaons or pions, up to two neu- 



K+TT-TT-TT^, K°Tr-, K^TT-TT^, K°Tr'TT~Tr+. Wc define 
the purity of a particular mode as S/{S + B), where S 
(B) denotes the number of signal (background) events; 
we use only the 186 -Broco final states with purity, mea- 
sured in data control samples, greater than 0.2. When 
more than one i?rcco candidate is found in an event, the 
one with the highest purity is retained; the overall purity 
of our selected sample is approximately 0.45. 

Two kinematic variables characterize correctly recon- 
structed B candidates: the energy-substituted mass 
Toes = y/s/4: — and the energy difference AE = Eb — 
-v/s/2, where {Eb,Pb) is the B-meson four-momentum 
in the T{AS) rest frame. For the Bj-eco candidate, we 
select events with 5.2500 < toes < 5.2893 GeV and we 
apply a mode-dependent cut on AE'. Additional back- 
ground rejection is provided by the angle 6t, defined as 
the angle between the thrust axis of the i?roco candidate 
decay products and the rest of the event. For contin- 
uum events | cos6't| peaks sharply at 1, while BB events 
exhibit a uniform distribution. We select events with 
|cos6't| < 0.9. 

Finally, we combine into a Fisher discriminant T four 
variables sensitive to the event shape and the production 
dynamics: the polar angles with respect to the beam axis 
in the T{AS) frame of the -Brcco candidate momentum 
and of the i?rcco thrust axis, and the zeroth and second 
angular moments L0.2 of the energy flow. The moments 
are deflned by Lj = J^iPi ^ |cos6'i|'' , where i labels a 
charged or neutral candidate not originating from the 
decay of the B^eco, Oi is the angle with respect to the 
-Breco thrust axis, and pi is its momentum. 

The branching fractions we are measuring are nor- 
malized to the number of fully reconstructed BB events 
present in our sample. We determine the BB yield (over 
the qq continuum background) through a maximum like- 
lihood fit to the variables toes and T. The probability 
density function (PDF) of toes for the BB category is 
the sum of two components: two Gaussian functions cen- 
tered on the mass of the B parameterize the correctly re- 
constructed B candidates, while an ARGUS [13] function 
describes the misreconstructed B decays. For the contin- 
uum we use only an ARGUS function. For the J- variable 
we use the sum of a bifurcated Gaussian with a Gaussian 
for both BB and qq. Besides the yields of the two com- 
ponents {BB and qq), the ARGUS exponent for the qq 
component and the fraction of correctly reconstructed 
BB events are free. We split the data sample into four 
subsamples characterized by different purity ranges of the 
Sroco candidates. The ARGUS exponent and the fraction 
of BB events peaking in toes are allowed to take different 
values among these categories. Figure [1] shows the pro- 
jection over the toes variable of this fit. The BB yield is 
(2.0902 ± 0.0020) x 10^ BB events. By repeating the fit 
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on the subsamples with different purities and using dif- 
ferent parameterizations for the PDFs, we estimate the 
systematic uncertainty on the BB yield to be 5%. 



xlO^ 




5.25 5.26 5.27 5.28 5.29 



iHes (GeV) 

FIG. 1: Projection of the tties variable for the Brcco sample; 
the dashed line represents the BB component, the dot-dashed 
is the continuum background, and the solid line is the sum of 
the two components. 



We assign to Bgig all the charged and neutral parti- 
cles not belonging to the Broco candidate and require 
5.1000 < mEs(-B,ig) < 5.2893 GeV. This loose cut sup- 
presses background events in which a significant amount 
of energy and momentum is lost. We suppress 6 — > c 
semileptonic decays by rejecting events where an elec- 
tron or muon candidate is present. We also veto events 
in which a D'^ , , or Df candidate, with a mass within 
30 MeV of the nominal value, is found. 

We require that a K^, K^, or 7r+ candidate with 
p* > 1.8 GeV be present on the signal side. The distance 
of closest approach for K'^ and tt"*" candidates must be 
less than three standard deviations from the Bgig decay 
vertex. K'^ candidates are reconstructed in the tt+tt^ fi- 
nal state, with requirements that the vertex probability 
of the two tracks be greater than 10~^, that the flight 
length be greater than three times its uncertainty, and 
that their mass satisfy 0.486 < m^+^- < 0.510 GeV. 

We extract the signal yields from a maximum likeli- 
hood fit to the three variables jties (-Broco), ^, and p* . 
For the and Tr~^ samples we also measure the direct 
CP asymmetry Ah = (r^ - F+)/(r- F+), where the 
superscript to the decay width F refers to the charge of 
the light hadron. Our fits have three components: sig- 
nal, 6 — >■ c background, and continuum background. For 
each of these categories j we define probability density 
fimctions Vj{x) for the variable x, with the resulting like- 
lihood: 



V, - V,{m^s)V,{T)V,[p*), (1) 

_ y. N 

^ - ^^^E^.n^ (2) 

1=1 j 

where Vj is Vj evaluated for event i, Yj is the yield for 
category j, and N is the number of events entering the 
fit. We assume the PDFs for each variable to be uncor- 
related in the signal and 6 — ^ c components (a correla- 
tion in the continuum component is handled as discussed 
below). We check this assumption by means of Monte 
Carlo (MC) experiments |l^, in which signal and 6 — c 
events are taken from fully simulated event samples and 
the continuum background is generated from the PDFs. 
In the extraction of the signal yields, we correct for the 
small biases we observe in these ensembles. The PDFs 
are extracted by fitting MC samples, where the charm- 
less decays are separated from 6 — > c background using 
information at the generator level. 

Signal and 6 — c events share the same PDFs for the 
Toes and variables which are only effective to separate 
BB events from the continuum; the fit distinguishes be- 
tween charmed and charmless B decays by exploiting the 
differences in the p* distributions. The p* distribution is 
parameterized by the sum of a Gaussian with an ARGUS 
component for the signal, by the sum of an exponential 
and a Gaussian for the qq component, and by the sum of 
three, one, or five Gaussians for the b ^ c background in 
the and 7r+ samples, respectively. The latter 

parameterize the broad component (s) of the b c back- 
ground and the peaking components corresponding to the 
B £)(*.**) /i, (h = K+,K° or 7r+) decays, ah of which 
are evident in the 7r+ sample (see Fig. [5]). Similarly, 
the Gaussian component of the signal p* PDF accounts 
for the dominant two-body decays (mainly B — )■ rj'K), 
while the broad component describes the sum of the other 
contributions. The splitting of the data into subsamples 
based on the purity and the charge of the -Brcco candi- 
dates allows differences in the background distributions 
to be accommodated in the fit by allowing the param- 
eters most sensitive to these variations to take different 
values in each subsample. 

The fit is performed through an iterative procedure. In 
the first step we fix the signal yield to the predictions of 
the MC and fit the p* > 1.8 GeV sample, leaving free to 
vary the most important parameters of the background 
such as the normalization of the peaking components in 
the 6 — ?> c background, the width of the broad compo- 
nents, and the exponent of the ARGUS function. This 
step is aimed at determining the shape and the normal- 
ization of the 6 — >■ c background; the projection plots for 
this step of the fit are presented in Fig. [2] 

In the next step, we use the results obtained in the pre- 
vious fit to extrapolate the predicted 6 — c background 
into the high p* region {p* > 2.34 GeV for and K^, 
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FIG. 2: Projection plots for the whole p* range for the (a) K'^ , (b) Kg, and (c) tt^ samples. The solid curves are the total fit 
functions, the red dashed lines are the signal components (which are kept fixed at this stage), the blue long dashed lines are 
the fe — >■ c background and the magenta dotted lines are qq. The scale on the upper border of the plots indicates the mass of 
the system recoiling against the light hadron. 



TABLE L Summary of the fit results to the high p* range. The b ^ c background yield is kept fixed in this fit; the quoted 
uncertainty represents the amount by which this quantity is varied for the evaluation of systematic uncertainties. The first 
error in the branching fractions and in the direct charge asymmetries is the statistical one, while the second is systematic 
(the significance includes only the additive part of the latter). The upper limits (U.L.) on the partial branching fractions are 
taken at the 90% confidence level. For the tv^ sample, the results of the yields refer to the p* > 2.36 GeV range, whereas the 
branching fraction has been extrapolated to p* > 2.34 GeV. 



B K+X 



B K"X 



B TT + X 



Events to fit 
6 — ^ c yield (events) 
qq yield (events) 
Signal yield (events) 
Fit bias (events) 
Significance (a) 

B (Xl0~ )p->2.34 GcV 

B U.L. (xlO-^ 



>2.34 GcV 



306 
66 ±8 
194 ± 15 

+10.9 
2.9 

"1+32 

-29 

187 



84 
6.5 ±2.6 

48 ±8 

+3.5 
3.8 

-+51 

-45 

294 



692 
173 ± 13 
430 ± 22 

-4.3 
6.7 

1+50 



0.57 + 0.24 + 0.05 



0.10 + 0.16 + 0.05 



p* > 2.36 GeV for 7r+). We fit these subsamples, varying 
only the yields of the signal and qq background compo- 
nents and the charge asymmetries, while the shapes are 
those determined in the previous step (see Fig. [3]). An 
exception occurs for the T variable in the qq background, 
which is correlated with p*; thus, fixing its shape to that 
determined in the whole p* range would lead to a bias. 
In this case we parameterize the jF distribution with two 
Gaussians, determine its parameters from the MC in the 
high p* range, and leave the mean of the core Gaussian 
free to vary in the fit. Using the p* cut efficiency de- 
rived from the MC, we then recalculate the number of 
signal events in the whole p* range and repeat the fitting 
procedure from the beginning. 

We find that this procedure converges after at most six 
cycles and that the result does not depend on the initial 
values we choose for the signal yield. We use the results 



of the final fit to the high p* range to derive the partial 
branching fractions and the direct CP-asymmetries (for 
the and 7r+ samples). The branching fractions are 
computed using the efficiencies for reconstructing signal 
events in the high p* region derived from the simulation. 
In order to avoid the systematic uncertainty related to 
the -Brcco reconstruction efficiency, the calculation is done 
taking for the normalization the number of BB events 
present in our sample. To make the comparison with the 
kaon samples easier, we extrapolate the branching frac- 
tion of B — >■ TT+X to the p* > 2.34 GeV range (we assume 
the systematic error associated with this extrapolation to 
be negligible) . The results are collected in Table HI 

The whole fit procedure is tested on a data sample en- 
riched in 6 — >■ c background, selected by reversing the 
vetoes on the 1?°, D^, or Df candidates associated with 
the i?sig- The results agree within statistical uncertain- 
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FIG. 3: Projection plots for p* > 2.34(2.36) GeV for the (a) K'^ , (b) K^, and (c) tt^ samples. The solid curves are the total 
fit function, the red dashed fines are the signal component, the blue long dashed are the 6 — >■ c background and the magenta 
dotted are qq. In order to enhance the signal component we apply cuts on the likelihood (computed excluding the p* variable) 
which retain 82 -88% of signal events while suppressing most of the qq background. The scale on the upper border of the plots 
indicates the mass of the system recoiling against the light hadron. 



ties with the expectations of very small signal yields. We 
also verify that our model for the continuum background 
is in very good agreement with the data taken away from 
the T{AS) resonance. 

Systematic uncertainties arise from the imperfect 
knowledge of the number of Broco candidates (5%), from 
the uncertainties on the reconstruction efficiencies for 
charged particles (0.5%), K'^ candidates (2.1%), and 
other neutral particles (0.9-1.2%, depending on the fi- 
nal state), from the K/tt separation (2.4%), and from 
the statistics of the MC sample which we use to compute 
the efhciency in reconstructing signal events (6.8-14.5%). 
The above uncertainties are multiplicative and do not af- 
fect the significance of the measured branching fractions, 
contrary to the following additive contributions: the un- 
certainty on the PDFs of the signal component is esti- 
mated by leaving each parameter kept fixed in the nom- 
inal fit free to vary (3.6-8.5 events). The uncertainty in 
the 6 — > c background is computed by varying its yield by 
the sum in quadrature of its Poisson uncertainty and the 
uncertainty in the extrapolation to the high p* region, 
taking into account the uncertainty on the knowledge of 
the signal PDF. The resulting systematic error is 2.8-10.3 
events. The systematic error arising from the correction 
for the fit bias is taken as the sum in quadrature of half 
the correction itself and the statistical uncertainty on the 
correction (3.6-7.9 events). 

The systematic uncertainties for the direct CP- 
asymmetries include the uncertainty in detector related 
charge asymmetries, which mainly affect the kaons (2%), 
different reconstruction efficiencies for B and B candi- 
dates in the tag sample (2.5%), and effects due to mistag- 
ging (3%). 

Our results for the partial branching fractions and Ach 



are given with statistical and systematic errors in Ta- 
ble m The central values for the branching fractions are 
in agreement with our estimates [3l of the sums of the 
known exclusive branching fractions of charmless two- 
and three-body i?-decays. On the other hand, predic- 
tions based on SCET underestimate the measure- 
ments, both those of the inclusive branching fractions 
presented here and those obtained by summing exclusive 
modes, even after adjusting for the branching fractions 
of the B — > -q^'^X modes, which are acknowledged to be 
problematic for SCET. This fact is interpreted by the au- 
thors of Ref. as an indication of the need to introduce 
substantial nonperturbative charming penguin contribu- 
tions or large higher-order corrections. 

In conclusion we have measured the inclusive partial 
branching fractions for B — >■ K^X, B — > K'^X, and 
B — > TT+X in the region where the momentum of the 
candidate hadron is greater than 2.34 GeV. The statisti- 
cal significance, computed as the difference between the 
value of — 21n£ for the zero signal hypothesis and the 
value at its minimum, exceeds five standard deviations in 
each case; however, comparable systematic uncertainties 
lower the significance to the values quoted in the Table, 
and we quote 90% confidence level upper limit (taken as 
the value below which lies 90% of the total of the likeli- 
hood integral, in the region where the branching fraction 
is positive) for the and K'^ modes. All results are in 
agreement with the standard model predictions, and ex- 
clude large enhancements due to sources of new physics. 
We do not find any significant direct (7P-asymmetry in 
the and 7r+ samples. 
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